Anomalous Density-of-States Fluctuations in Two-Dimensional Clean Metals 
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It is shown that density-of-states fluctuations, which can be interpreted as the order-parameter 
susceptibility xop m a Fermi liquid, are anomalously strong as a result of the existence of Goldstone 
modes and associated strong fluctuations. In a 2-d system with a long-range Coulomb interaction, a 
suitably defined xop diverges as 1/T 2 as a function of temperature in the limit of small wavenumber 
and frequency. In contrast, standard statistics suggest xop = 0(T), a discrepancy of three powers 
of T. The reasons behind this surprising prediction, as well as ways to observe it, are discussed. 

PACS numbers: 



It is known that the Fermi-liquid (FL) state can be de- 
scribed in terms of a stable renormalization-group (RG) 
fixed point (FP) Q. A related concept is the possibility 
of understanding the FL as an ordered state, in anal- 
ogy to the ferromagnetically ordered state and its cor- 
responding stable RG FP @, This notion was pio- 
neered by Wegner in the context of disordered systems 
More recently it was realized that the principle be- 
hind it is much more general and powerful, and can be 
applied to clean FLs as well 0, Important features 
of the resulting effectyive field theory are that it inte- 
grates out all massive degrees of freedom to arrive at an 
effective soft-mode theory that allows for a very straight- 
forward analysis of the FL FP, associated corrections to 
scaling, and the scaling behavior of various observables. 
Within this conceptual framework, the order parameter 
(OP) is the density of states (DOS) Q, the field con- 
jugate to the OP is a frequency-dependent chemical po- 
tential, the spontaneously broken continuous symmetry 
is a rotational symmetry in frequency space that can be 
understood as the symmetry between retarded and ad- 
vanced degrees of freedom, and the resulting Goldstone 
modes are the soft particle-hole excitations. These quan- 
tities are analogous to the magnetization m, the magnetic 
field h, the spin-rotational symmetry, and the magnons, 
respectively, in a Heisenberg ferromagnet. Note, how- 
ever, that the Goldstone modes in the FL case are soft 
only at zero temperature (T = 0). 

Let us recall the behavior of the DOS as a func- 
tion of the frequency uj (or the energy distance from 
the Fermi surface). For our purposes we mainly fo- 
cus on two-dimensional (2-d) systems with a long-ranged 
Coulomb interaction, which is arguably the most inter- 
esting case. However, the theory can be generalized to 
arbitrary d > 1, and the case of a short-ranged interac- 
tion is also of interest, see below and Ref. H|. It is known 
that the DOS is a nonanalytic function of uj, namely, 



a is a dimensionless coefficient. o(u>) denotes terms that 
vanish faster than linearly for uj — > 0. This result fol- 
lows from the soft-mode theory developed in Refs. 
(see below); technically, the frequency correction in Eq. 
([1]) is the leading correction to scaling at the FL FP. It 
has been calculated explicitly within many-body pertur- 
bation theory, which yields a weak-coupling value for a, 
viz., a = 1/4 It is important to note, however, 

that perturbation theory alone can never establish the 
\uj\ as the true leading frequency dependence. Our RG 
analysis of the soft-mode field theory, on the other hand, 
guarantees that the exponent of unity is exact. The pref- 
actor a will acquire corrections if one goes to higher-loop 
order. The linear frequency dependence is also consis- 
tent with simple scaling considerations as it should 
be It has been observed in high- mobility tunnel 



junctions [14 1 



N{uj) = 2AT F (1 + a \uj\/e F ) + o(u>) (1) 
Here Np is the DOS per spin at the Fermi level ep, and 



With these results as a backdrop, it is very interest- 
ing to ask about the behavior of the OP suceptibility. 
This quantity is an important measure of the strength 
of fluctuations in the ordered phase, which in turn pro- 
vides a measure of the likelihood that the ordered phase 
will become unstable as a function of a suitable param- 
eter. In the case of the classical Heisenberg ferromagnet 
mentioned above, the OP susceptibility is the magnetic 
susceptibility Xm = dm/dh. As a function of the wave 
number k, the longitudinal component of Xm diverges for 
k — > for all dimensions d < 4: Xm,h(k) oc l/k i ~ d [l5j |. 
This reflects the strong fluctuations in the ordered phase 
due to the existence of Goldstone modes that couple to 
the OP. It is then natural to ask whether an analogous 
susceptibility exists in the FL, and what its behavior at 
small wave numbers and frequencies is. 

To answer this question, we have performed explicit 
calculations [8|, some of which are quoted below. How- 
ever, for the sake of physical clarity it is very illustrative 
to derive the results from qualitative arguments, which is 
the purpose of this Letter. To this end, we first represent 
the DOS as 



N(u) = Re G(a; = 0, iu n -> uj + iO) /w 



(2) 
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Here 

G(x = 0,iuj n ) = (p(x,iu> n )) 
is the exact Green's function, with 

p{x, IbJn) = ^n(aj) 1pn(x) 



(3a) 



(3b) 



in terms of fermion fields t/) n (x) and ipn( x ) that depend 
on a fermionic Matsubara frequency u n = 2irT(n + 1/2), 
(n = 0, ±1, . . .). We suppress the spin degree of freedom, 
which is not important for our purposes. p{x, iuj n ) can 
be interpreted as the density of particles at point x and 
(complex) energy iuj n [HI- The temporal Fourier trans- 
form is defined as 



il> n (x) = T 



- T 1 / 2 



1/T 



dT 



tp(x,T) 



(4) 



with r the imaginary time variable. Dimensionally, 
ij)(x,T) is an inverse square root of a volume, so N is 
dimensionally indeed a DOS. Now we consider the spa- 
tial Fourier transform p(k,iu) n ) — J dx e~ lk ' x p(x,iuj n ), 
and the fluctuation Sp = p — {p}, and define a correlation 
function 



Xop(k,iu} n ) = (dp(k,iui n ) 5p(-k,iu n )) 



(5) 



Xop is the order-parameter susceptibility of interest. It 
describes both the response of the OP (i.e., the DOS) 
to it conjugate field, and the spontaneous fluctuations 
of the OP. The former is related to the latter by the 
fluctuation-dissipation theorem. 

Let us first use simple statistical arguments to deter- 
mine the behavior of xop one would expect in the absence 
of anomalous fluctuations. This will also serve to make 
a general point that is related to the definition of the 
temporal Fourier transform above. Consider <p n {x) = 
tp n (x)/VT, and define the "volume" V T = 1/T in the 
imaginary-time direction of the space-time of quantum 
statistical mechanics. Then (p n (x)) = (<p n (x) ip n (x)) oc 
Vt is a "time-extensive" quantity that is proportional 
to Vt = 1/T. Now consider the relative fluctuation 
((6p n (x)) 2 )/(p n (x)) 2 ex V T /V 2 = 1/V T = T. This just 
says that the relative fluctuation is proportional to 1/Vt, 
as one would expect from ordinary statistics. This yields 
an estimate for the fluctuations of p: 



((Sp n {x,iuj n )) 2 ) oc 



V 2 



({5p n ) 2 )/V 2 

oc 1 / Vt = 2 • 



(6) 



{(pn)/v T y 



All of these simple arguments of course assume that 
there are no strong fluctuations in the system that in- 
validate the simple statistics. Independent of this as- 
sumption, however, they show that ((Sp n (x,iui n )) 2 ), and 
hence XOP; must carry a factor of T that is simply the 
inverse linear extension 1/Vt of the imaginary-time axis. 
We thus have 



In the absence of anomalous fluctuations, 8 will scale as 
the zeroth power of the wave number, the frequency, or 
the temperature; i.e., 6 ~ 1. 

It is easy to see that the latter conclusion cannot be 
correct in an actual FL. Consider the scaling arguments 
from Ref . [l7|, and first consider the case of a short-ranged 
interaction, in which case there is only one dynamical 
scaling exponent z = 1 . This reflects the fact that in the 
Goldstone modes that characterize the broken symmetry 
in the FL, the frequency scales as the wave number, uj ~ 
k. At the stable FL FP, and in d = 2, the scaling part of 
the free energy density / obeys the homogeneity law 



f(T, h) = b- {2+z > f(Tb z ,hb z 



(8) 



Here we have assigned scale dimensions [L] = —1 and 
[T] = z to factors of length and temperature (or energy 
frequency or inverse time), respectively, b is the RG 
length rescaling factor, h is the field conjugate to the 
OP, which has a scale dimension [h] = z since the DOS 
must scale as an inverse energy times an inverse volume. 
Frequency-dependent quantities will scale with respect 
to frequency the same as with respect to temperature. 
Differentiating once with respect to h, we obtain for the 
scaling part of the OP density, i.e., the DOS correction 
SN = N - 2A F , 



6N(T) = b- 1 SN(Tb) cx T 



(9) 



in agreement with Eq. ([T]). Differentiating again with 
respect to h yields 



XO p(T, k) = X op(T6, kb) = f x (T/k) 



(10) 



with f x a scaling function. If we combine this with the 
conclusion from Eq. ([7]). namely, that xop must be pro- 
portional to T for fundamental quantum statistical rea- 
sons, we conclude 



Xop(T, k) oc 



T/k 



for T -> at fixed k 



T/T = 0(1) for k -> at fixed T 



Xop(k,iu n ;T) = T9(k,iuj n ;T) 



(7) 



, (11) 

This result is not consistent with the naive statistical 
arguments given above: If the DOS were normally dis- 
tributed, we would have xop oc 1/Vt and 6 ~ 1. In- 
stead, we see that the quantity 9 in Eq. ([7]) scales as 
6 ~ 1/fc ~ 1/T, and xop ~ 1- This divergence of the 
relative DOS fluctuations reflects the strong fluctuations 
in the system that are a consequence of the existence of 
the Goldstone modes. 

In order to understand what happens in the case 
of a long-ranged Coulomb interaction, we consider the 
structure of the field theory developed in Refs. 0, @. 
Within this framework, the DOS is schematically given 
by a frequency-momentum integral over a vertex func- 
tion V(k,iu>) times the potential U(k,iu>), with k the 
wave number and iuj the imaginary frequency. The vertex 
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function has the structure V(k,iu) = v(uj/vFk)/k 2 , with 
v(x — > oo) oc 1/x 2 . U is a constant in the short-range 
case, and represents the dynamically screened Coulomb 
potential in the long-range case. In d = 2, the latter has 
the schematic structure 

U(k,iuj) = 1 ' 



^/l+U J /(»Ffe) 2 



k + K,(vpk) 2 /ui 2 



(12) 

Here k is the screening wave number, and the second 
expression is valid in the limit v?k <C lj. In the denomi- 
nator we recognize the plasmon excitation, with its char- 
acteristic u) ~ fc 1 / 2 scaling. We also see the well-known 
fact that screening breaks down at nonzero frequencies. 
The DOS is schematically given by an integral (l9j 

/•A />ufA -i 

N ~ dkk du -^v(uj/v F k)U{k,iuj) , (13) 

Jo Jt k 

with A an ultraviolet momentum cutoff. Power count- 
ing yields N(T) oc const. — T, in agreement with Eqs. 
(j9]) and ([T]). This is true in both the short-range and 
the long-range case; the singular nature of the screened 
Coulomb interaction at nonzero frequency does not suf- 
fice to change the behavior in d = 2 20] . 

Now consider the OP susceptibility. It follows from 
Eq. §5§ that xop is given by an integral analogous to the 
one in Eq. fjl3|) with the integrand squared and an extra 
factor of T: 



Xop 



T 



dk k 



v F A 



du> ( —v(ui/vFk)U(k,iuj) 



(14) 

The prefactor of T is just the trivial "volume" factor 
1/Vr discussed in the context of Eqs. ((HIT]). In the short- 
range case, we obtain xop oc T/T in agreement with Eq. 
(fTTj) . In the long-range case, a power-counting analysis 
of the integral yields xop °c T/T 3 , with a multiplicative 
logarithmic infinity coming from the momentum integral. 

A careful analysis of xop within the field theory of 
Ref. 0, or its generalization to a long-ranged Coulomb 
interaction, yields results that are consistent with the 
ones sketched above. For the static OP susceptibility in 
the limit of low temperatures we find, in the long-range 
case, 

2 rp 

XOP(fc) = ^br T 3 ( 15a ) 

This result is valid for T 2 /vpk 3> v-pK, and to logarithmic 
accuracy. In the opposite limit, we find 

k 1 / 2 T 

Xop(k) = const, x —3- , (15b) 

with const, a constant of O(l) [2l[. Finally, for T — >• 
and nonzero external frequency ui one finds, to logarith- 
mic accuracy for uj 2 3> v^nk, 

8 n 2 T 

Xop(k,iuj) = —r — 3- ln(w 2 /u|Kfc) (15c) 
6ir z oj a 



As expected, the frequency scaling is the same as the 
temperature scaling. Comparing Eqs. (115b[) and (|15a[ 
I15cj) we see that Toru scale as A: 1 / 2 , which reflects the 
fact that the integral is dominated by the plasmon time 
scale. This is also apparent from the arguments of the 
logarithms in Eqs. (|15aU15cj) . The fact that xop(k — > 0) 
diverges even at a nonzero T is due to the fact that the 
plasmon, being a density fluctuation governed by a con- 
servation law, remains soft even at T > 0. Consequently, 
there is a logarithmic divergence (in d = 2) that is not 
cut off by temperature. 

Also of interest is the homogeneous susceptibility in a 
finite system of linear dimension L. In that case, k in 
Eqs. (fl~5|) gets replaced by 1/L, and we have 



for T 2 L > u 2 k 



Ul/T 2 )\n{T 2 L/v 2 K) 
XOP oc I TL^K l ' 2 /vl 

((T/uj 3 )Itl{u 2 L/v 2 K ) 



for T 2 L <C vIk 

for T -t 0, uj 2 L > v%k 



(16) 

In summary, we have calculated the susceptibility that 
describes density-of-states fluctuations in a 2-d Fermi liq- 
uid; a correlation function that can be interpreted as the 
OP susceptibility of the ordered phase represented by the 
Fermi liquid. The OP susceptibility, properly normalized 
to account for a trivial factor of T, see Eq. (JT]), diverges 
in the limit of vanishing wave number k as \/k 3 / 2 , Eq. 
(|15b|) . or as 1/T 3 or 1/lj 3 in the limits of vanishing tem- 
perature or frequency, respectively. This is a very strong 
effect; uncorrelated statistics would lead to a constant 
for the same quantity. It is the combined result of the 
Goldstone modes, i.e., particle-hole excitations, in the 
ordered phase, and the plasmon excitations in an elec- 
tron system with long-range Coulomb interaction. The 
frequency scale associated with the latter (u> ~ fc 1 ^ 2 ), 
although formally irrelevant compared to the Goldstone 
mode frequency scale (u ~ k), acts as a dangerous ir- 
relevant variable in the RG sense and leads to an even 
stronger divergence than one would expect from naive 
scaling. Even in the case of a short-range interaction 22j, 
where only the Goldstone modes are present and naive 
scaling works, the normalized susceptibility diverges as 
1/k ~ 1 jT ~ 1/w, compared to a constant in the absence 
of anomalous fluctuations. A simple analogy of this ef- 
fect in a classical system is the l/k A ~ d divergence of the 
longitudinal spin susceptibility in the ordered phase of a 
Heisenberg ferromagnet. 

Physically, this anomaly is very closely related to the 
|o;| anomaly in the DOS, Eq. ([T]), which constitutes a 
pseudogap as originally defined by Mott [23| ■ It is worth 
noting that this pseudogap is an intrinsic feature of any 
2-d Fermi liquid. We stress again that the current theory 
shows the linear frequency dependence to be exact, which 
cannot be established within perturbation theory. 

We finally add some comments about the experimental 
relevance of the above results. The spectrum of the local 
Green's function, Eq. (|3"a|) . is what is often referred to as 
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the local density of states (LDOS) [24j,l25|. The LDOS 
gives the dominant contribution to the tunneling current 
in a scanning tunneling experiment [26]; its spatial av- 
erage, which is the DOS, is measured in a classic tunnel 
junction. Our OP susceptibility, Eq. ([5]), describes the 
spatially averaged fluctuations of the LDOS. A two-tip 
tunneling experiment 0, H?} should be able to give in- 
formation about this quantity jjjj ]. 
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